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Abstract 



The large time and long distance behavior of the temperature correlation functions of the 
quantum one-dimensional Bose gas is considered. We obtain integral equations, which solu- 
tions describe the asymptotics. These equations are closely related to the thermodynamic 
Bethe Ansatz equations. In the low temperature limit the solutions of these equations are 
given in terms of observables of the model. 

1 Introduction 

This paper continues the series of works 0-fl, devoted to the study of the correlation functions 
of the quantum one-dimensional Bose gas with delta-function interaction. Here we consider the 
problems related to the large time and long distance asymptotics of the correlation functions. 
The main result of this paper is a system of integral equations (see section 4), which solutions 
describe the asymptotics. The form of these equations is close to the thermodynamic Bethe Ansatz 
equations, and apparently their solutions are closely related to observables of the one-dimensional 
Bose gas model. Anyhow such the relationship does exist in the low temperature limit. 

The main object of our investigation is a temperature correlation functions of local fields 



where c > is the coupling constant, h — chemical potential. The operators ty(x,t) and ^(x,t) 
are the canonical Bose fields 




tr e t 




Here T is a temperature, H is the Hamiltonian 
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For the reader's convenience we present here the list of the basic equations for observables of the 
one-dimensional Bose gas at a finite temperature. The energy of elementary excitation (particle) 
satisfies the Yang- Yang equation || 

e( A) = A>-k--/_X(A, < 0to(l + «--* 1 )4,, K(X,,) = ^— ? . (1.4) 

In the case of positive chemical potential the function e(X) has two real roots e(±gy) = 0. The 
total spectral density of the vacancies in the gas is given by the equation 

/oo 
K{\^)^)p t {^)d^ (1.5) 
-oo 

where 



0(A) = 1 + exp 



e(A) 
T 



;i.6) 



plays the role of the Fermi weight. 

The momentum of particle, being a function of the spectral parameter, also can be found 
from the corresponding integral equation. However below we shall need only the derivative of the 
momentum with respect to the spectral parameter. This quantity coincides up to a coefficient 
with the total density of the vacancies: 

^ = 2 W (A). (1.7) 

The velocity of the particle v(X) is 

, , de e'(\) e'(A) , , 

Finally, we present also the integral equation for the scattering phase 

/OO / jp _|_ \ — y\ 
K(\u)$(u)F(u,v)du + i\n[ — . (1.9) 

-oo \ic + v — A J 

Here and further the branch of a logarithm is chosen such that — 7r < arglnz < tt. The quantity 
5(A, v) = exp{27riF(A, u)} is the scattering matrix of two particles possessing momenta fc(A) and 
k(v). 

All the thermodynamics observables listed here will be needed later on for the description of 
the correlation function asymptotics. 

We shall consider the behavior of the correlation function ( p..l|) at x — ► oo, t — > oo and fixed 
ratio Ao = x/2t. The one- dimensional Bose gas is massless model, therefore the correlation function 
( |1 . 1| ) decay as a power at zero temperature 

(^(o,o)^(x,t)) ^cr Ao . 
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The exponent A can be computed in the framework of the conformal field theory J7[|. At finite 
temperature the asymptotics replaces by the exponential one 



(*(0, 0)^(x, t)) T — ► Ct~ A e~ t/r . (1.10) 

In other words the zero temperature turns out to be the point of phase transition. The methods of 
the conformal field theory can be applied at a finite temperature only approximately, for instance, 
for calculations of the low temperature contributions into the correlation radius [0. 

The determinant representation method, used in allows one in principle to consider 

correlation functions at arbitrary temperature. This method is especially powerful for study of 
their asymptotics. In the present paper we focus our attention at the main, exponential law of 
decay in the equation (|1.10|) . We shall not concern neither the pre-exponential factor t~ A , nor the 
constant coefficient C in ( |1.10|) . In other words our goal is to compute the correlation radius. 



Let us give now the general content of the paper. In the next section we present a brief 
description of the method, based on the representation of the correlation functions in terms of 
Fredholm determinants. Here we also discuss the questions related to the construction of the dual 
fields. In particular we explain the cause of these auxiliary quantum operators appearance and 
their role in the calculation of the correlation functions. On the whole the section 2 is a sort of 
review. The reader can find the details in the papers mentioned above. In the section 3 

the methods of the averaging of the operators, depending on dual fields, is developed. Here we 
define the 'averaging mapping', associating quantum operators with certain classical functions. 
This mapping allows us to compute the asymptotics of the correlation function (|1 . 1|) . This is done 
in the section 4. Here the system of integral equations, describing the asymptotics, is given. In 
the last section we consider some specific cases of the correlation function (|1.1|) . In particular we 
prove that in the low temperature limit our results exactly coincide with the results of the paper 



2 Dual fields 

The method of the Fredholm determinant representation for evaluation of the correlation functions 
is described in details, for example, in the book ||. At the first stage of this method one need 
to represent a correlation function in terms of a determinant of an integral operator, which kernel 
depends on distance and time (and other physical parameters of the model). For the correlation 
function (|1 . 1| ) such the representation was found in jl]]. Further investigation of the Fredholm 
determinant obtained, can be performed via the methods of classical exactly solvable equations. 
The matter is that the determinants mentioned above turns out to be r-functions of classical 
integrable systems. In particular, for the model of the quantum one-dimensional Bose gas they 
can be expressed in terms of solutions of classical forms NLS (in the simplest case of two-points 
field correlation function in the free fermionic limit the corresponding differential equation is 
scalar Nonlinear Schrodinger equation). In turn the solutions of classical integrable equations can 
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be found by means of the Riemann-Hilbert problem approach. In particular, for the calculation 
of the large time and long distance asymptotics the nonlinear steepest descent method is used 
0, allowing in principle to obtain the complete asymptotic expansion for correlation functions. 
Thus, the computation of the correlation function asymptotics in the framework of the method 
described, consists of two stages: 

1) representation of the correlation function in the Fredholm determinant form; 

2) computation of the asymptotics of the determinant obtained via the methods of the Rie- 
mann-Hilbert problem and classical exactly solvable equations. 

However that is how the matter stands only for the free fermionic models (in the quantum NLS 
this corresponds to the limit of infinite coupling constant: c — ► oo). Out off the free fermionic point 
the situation is more complicated, and one more step should be added to the scheme described 
above: 

3) averaging with respect to an auxiliary vacuum. 

The present paper is devoted just to the last stage of calculations, therefore we consider this 
problem more detailed. 

In order to obtain the Fredholm determinant representation for the correlation functions of 
non-free fermionic models (finite coupling constant for one-dimensional Bose gas), one need to 
introduce auxiliary quantum operators — Korepin's dual fields [fLO) . This necessity is caused by the 
existence of a non-trivial ^-matrix. Consider the main idea in brief. For evaluation of a correlation 
function of an operator O with respect to a state \a) one can use the form factor decomposition 

(a\O(x,t)O(0,0)\a) = J2\(a\O(0,0)\b)\ 2 e D ^. (2.1) 

Here the state \b) runs through a complete set, the function D a b(x,t) describes the dispersion law 
of the model. The quantity (a\O(0,0)\b) is refered as a form factor of the operator O, and, for 
example, in the quantum NLS it can be computed explicitly via the algebraic Bethe Ansatz. Thus, 
for O = \I/ the form factor is proportional to a certain determinant: 

(a|tf(0,0)|6) ~ detM. (2.2) 

We would like to emphasize that in distinction of the quantum field theory models here the form 
factors are constructed in the bare Fock vacuum, but not in the physical one. If the states \a) and 
\b) are parameterized by the spectral parameters {A} and {/i} respectively, then the entries Mj k 
are functions of these spectral parameters 

M jk = M jk ({X},{fi}). 

The type of dependency of Mj k on the spectral parameters is extremely important. In fact, in 
order to obtain a determinant representation for the mean value ( |2.1| ), we need to reduce the 
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infinite sum of determinants in the r.h.s. ( |2.1| ) to a single determinant. This can be done only for 
very special form of matrices M, namely, if Mj k = M(Xj,[ik), i. e. the entries are parameterized 
by a single two- variable function M(x,y), taken in the points x = Xj and y = These are so 
called local matrices. In the reality the entries of the matrix Mj k , describing the form factor fl2.2|) , 
contain the products of two-particle S'-matrices 

S(A,„) = -*i±Az£ (2.3) 

with respect to all parameters {A} and {/i}. This leads to the effect that for interacting fermions 
the entries Mjk depend on all spectral parameters Mj k = M({A}, {/i}) (non-local matrices). As 
a result, it becomes impossible to reduce the sum of determinants ( |2.1| ) to a single determinant. 
On the other hand, in the free fermionic limit (c — > oo) the S-matrix becomes trivial, and, as a 
consequence, the entries Mj k depend only on two variables Mj k = M(Xj,fi k ). The details and all 
explicit equations can be found in ffl]. 

The introducing of dual fields updates the situation. In order to reduce the matrix M to the 
local form it is sufficient to factorize the S'-matrix 

S(A,/i) = A(A)/ 2 (/i), (2.4) 



where fj are some functions. Of course, such a factorization is impossible for the S"-matrix fl2.3|) , 
if fj are classical functions. However one can easily achieve this by means of quantum objects. 

Consider two creation operators q^p(X) and %,(A), and two annihilation operators p^(X), p</>(A), 
acting in an auxiliary Fock space as 

((%(A) = (O|g (A) = 0, ^(A)|0) = p (A)|O) = 0. (2.5) 

The non-vanishing commutation relations are 

MA), QM] = -M A ), qM] = In ( ^^-^j j ■ (2-6) 

The dual fields are linear combinations 

V>(A) = <fy(A) +^(A), 0(A) = q^X) +p^(X). (2.7) 

Obviously, due to ( |2.6|) the dual fields commute with each other 

[V(A), il>(/*)] = [0(A), 00*)] = [V(A), 0^)] = 0- (2.8) 

On the other hand, the vacuum expectation value of the expressions, containing the dual fields, 
might be non-trivial 

S (X,fi) = ^ C + X ~f i = _( |e*W e ^)|o). (2.9) 
ic + n — X 
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Thus, the S'-matrix can be presented as vacuum expectation value of the factorized expression of 
the form (|2.4j) . Due to the property ( |2.8|) , the determinants, depending on the dual fields, are well 
defined. As a result the form factors can be presented in the form 

<a|O(0, 0)|6> ~detM= (0| detM(0, ip)\0), (2.10) 

where matrix M is now local one, i. e. Mjk = M(\j,[ik)- 

Similar representation can be obtained for the square of the absolute value of the form factor 
(|2.1|) , and the summation with respect to complete set {\b)} can be performed under the symbol 
of the vacuum expectation value. Eventually we arrive at the representation for the mean value 
( |2.1| ) in terms of the vacuum expectation value in the auxiliary Fock space of a determinant of a 



certain matrix, depending on the dual fields. In the thermodynamic limit this matrix turns into 
an integral operator, and we obtain 

(C(x,t)C(0,0)) ~ (0|det(/ + y)|0). (2.11) 



Here the kernel of the operator V depends on the dual fields 0(A) and ^(A), however due to 
the Fredholm determinant is well defined object. 

We would like to emphasize that here we have described only the main idea of the dual fields 
method. In the specific cases some modifications are possible, in particular, sometimes one need 
to introduce additional dual fields. The reader can find the detailed formula? in jl]. However the 
basic idea of the method always consists of the applying of factorization (|2.4|), which allows to 
reduce the determinants of non-local matrices to the local ones. 

Thus the correlation function of the quantum one- dimensional Bose gas turns out to be pro- 
portional to the vacuum expectation value of the Fredholm determinant, functionally depending 
on the dual fields. Nevertheless the asymptotic analysis of this determinant can be performed in 
the same manner, as for free fermionic limit. Due to the commutation relations ( |2.8|) at the certain 
stage of the calculations one can consider the dual fields as some classical functions, which are 
analytical in the strip < c/2. The last property follows from the representation of the dual 
fields in terms of the canonical Bose fields 



1 r ( 1£ _|_ I) — \ \ 

^ y ! VlnJ (*/-A) 2 + f 



(2.12) 



where 

[<Pi(\),<pU»)] = [<P2WM{»)]=8(\-n), ^(A)|0) =0, (Oj^t(A) = 0. (2.13) 

It is easy to check, that commutation relations between the operators and q^, defined by 
( 2.12 ), exactly coincide with ( |2.6| ). 
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Via the Riemann-Hilbert problem and classical integrable equations methods one can find 
the large time and long distance behavior of the Fredholm determinant fl2.11|) (see ||-f5[]). The 
asymptotic expression obtained is a functional of the dual fields, which up to this moment are 
considered as some classical functional parameters. In order to calculate the asymptotics of the 
correlation function one need now to remember about quantum nature of these operators and to 
average the obtained expression with respect to auxiliary vacuum. 

Here, however, one important moment exists, which had been discussed in details in 0. The 
matter is, that the procedures of averaging and calculation of the asymptotics, strictly speaking, 
do not commute with each other. The asymptotics of the correlation function is equal to the 
asymptotics of the vacuum expectation value ( |2.11| ), but not to the vacuum expectation value of 



the asymptotics. Just because of this reason the results of the work []TT| , where some questions 
of averaging of the expressions containing dual fields, were considered, can be applied to the 
evaluation of the correlation functions only in some measure. Nevertheless it was proved in the 
mentioned above paper || that there exist asymptotic representations, which are stable with 
respect to the procedure of averaging, i. e. 

asymptotics(vacuum expectaion value) = vacuum expectation value(asymptotics). 

Such the representation for the Fredholm determinant was obtained in M, and further we shall 
deal just with this one. 



3 Averaging mapping 

The asymptotic formulae, obtained in || for the Fredholm determinant in terms of the dual fields, 
are rather bulky (see (|4.1| ), ( f4.2|) below). Therefore before averaging of these expressions we con- 
sider several more simple, but at the same time quite general examples. The technique, developed 
in the present section, will allow us to average the equations ( |4.1| ), ( |4.2j ) easily. 

Consider two dual fields <p and ip, defining by (|2.7|). Let the commutation relations between 
the creation and annihilation operators are 

where £(A,/z) is a two- variable function. Here we do not need the explicit form of this function, 
however it is easy to see that putting £(A, /i) = ln(zc + fi — A) — ln(zc + A — //), we reproduce the 
commutation relations (|2.6|) . Besides the relations ( |3.1|) we introduce an additional one 

[?ty(A), tyifi)] = v(\ V) = V(V, A), (3.2) 

where fi) is a symmetric two- variable function. Below we shall see that the additional relation 
( |3.2| ) does not affect essentially on the final result, however for some time we do not put the 
function 7/(A, /i) equal to zero. 
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is still valid. However the 



It is easy to see that the main property of the dual fields 
vacuum expectation values may be non-trivial as before. 

The main goal of this section is evaluation of the mean value of the following form: 

( O | e ^)F(0)|O). 



(3.3) 



Here J 7 (if), <fi) and F(<f>) are operator-valued functionals (later on functionals for brief) of the dual 
fields. Functional F does not depend on the field tf>, while T depends on this field linearly 



b 2 T 

5ip 2 



0. 



(3.4) 



All the functionals (and functions) of the dual fields should be understood in the sense of formal 
series. 



We are starting with the formula, obtained in [11] 



(O|e)q>(x;a^(A fc )/ fc (0^ 

U=l J Sn=l ' aet 9 



where 



N 



E{7j) = exp \ - Z nZmV(K, K 

n,m=l 



(3.5) 



(3.6) 



Here in the l.h.s. f(^)(v)^ is a functional (or a function) of the field <p, a is a complex parameter, 
functions fk(z) assumed to be holomorphic in a vicinity of the origin. In the r.h.s. the quantities 
Zj are the roots of the system 



N 

: ; - a fj( Yl Am •/'.,) ) = 0. 

m=l 



(3.7) 



If the system has several solutions, then one should choose the single one, approaching zero at 
a — > (see ||11|| ). Finally, in the denominator of (3.5) the Jacobian of the system ( |3.7| ) is placed 



9jk 



_d_ 



V m=l J 



(3.8) 



It is suitable to modify the result (|3.5| ). Namely, let us introduce a function w(u) 

N 

w(u) = z m ^(\ m ,u). (3.9) 

m=l 

Multiplying each of the equations ( |3.7| ) by £(\j,u) and summing up all of them, we obtain 

G(w) = w(u) - a Y fm(w(fi m ))^(Xm,u) = 0. (3.10) 



m=l 



Thus instead of the system ( p.7| ) we obtain only one equation ( 5TTD ) for the function w(u), which 
has to be valid, however, for any value of u. The roots of the system Zj can be expressed in terms 
of w(u) as follows 

Zj = afj(w(jij)). (3.11) 

It is easy to see also that 

~6G(w(u)Y 



det gjk = det 



Sw(u') 



(3.12) 



Thus the result (B.5I) turns into 



v 



exp \ £ a^\ k )f k {<P^ k )) \ F(<p(u))\0) = F (w(v))-^. 

U=l J aet Sw 



(3.13) 



Let now the functional F(<f>) is simultaneously a function of some parameter t: F((j>) = F(<f>\t). 
This parameter will play the role of time in the asymptotic formulae for the correlation function. 
Then the equation ( |3. 13|) can be written in the form 



exp [j2 ^(h)fk(4>(Vk))\ F(<f>{v)\t)\0) = CF(w(u)\t), 



(3.14) 



where C is a constant, does not depending on t. Constant factors of the such type will be not 
interesting for us in computation of the correlation function asymptotics. Just because of this 
reason on can put the function rj(\, fj,) equal to zero, since it contributes only into the common 
constant multiplier and does not affect on the dependency of the r.h.s. of ( |3.14j ) on t. 

Thus we see that under the averaging of the expressions of the form ( p.5| ) a mapping of the 
operator <f) into the classical function w takes place 



<j){v) -> w(v), 



v 



F(w(i 



(3.15) 



where w(y) is the solution of the equation ( p. 10 ). We call the mapping ( 3.15 ) (as well as more 
general equations ( |3.17|) , (|3.18|) ) averaging mapping. 
Consider now a generalization of the functional 



N 



a 



k=i 



(3.16) 



standing in the exponent of ( |3.14 ). It is clear, that one can easily include the derivatives of the 
fields ift( m \X) and 0^ m ^(A) into this functional. One can consider also the continuous limit, when 
the sum with respect to k turns into an integral (this was done in |T]J). Besides it is possible 
to consider infinite series instead of finite sum with respect to k. In all listed cases the result of 
averaging can be easily obtained by the corresponding limiting procedure in the equation (|3.14|) 
(or in (|3.13|) , if more precise results are needed). It is important only that T{^), (ft) would be always 
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linear with respect to the field if) and its derivatives, while the dependency of T and F on the field 
(f> might be rather complicated. 

Thus we arrive at the following formula 

(O| e WW)F(0(*/)|t)|O) = CF(w(u)\t), (3.17) 
where the function w(u) satisfies the following equation 

w(u)-J 7 ((-(;u),w) = 0. (3.18) 

As for functionals, depending on the field if) in non-linear way, such a simple formula, generally 
speaking, does not exist. However in scattered instances it is possible to generalize the averaging 
mapping for more complicated functionals. Below we consider an example of functionals, which 
are 'asymptotically linear' with respect to the field if)(\). 

Consider a mean value 

(O|e^(-)^ A )F c (0(z/),A)e* Fe ^ M ' A )|O), (3.19) 

We are interesting of the behavior of this mean value at t — > oo. Let the explicit dependency of 
the functional ^F(ifj(-),(f),A) on the field if> is linear as before. Besides let T and F depend on if) 
implicitly through an operator A, which is defined as the 'root' of the equation 

A = Ao + ^'(A), A^Ao at t -> oo. (3.20) 

Here Ao is a constant, prime means the derivative with respect to the argument. The definition 
(|3.20|) is formal, and one should understand it in the sense of a series with respect to t -1 . If we 
replace in ( |3.20| ) the operator if) by a complex function holomorphic in a vicinity of Ao, then for 
arbitrary function f(z), holomorphic in the same vicinity, we have under the corresponding choice 
of the contour 

1 I 



/(A) = — f dzf(z) ^ 77? v ( 3 - 21 ) 

2m J z - A - j- t tp [z) 

For large enough t one can expand the denominator in (|3.21| ) into the absolutely convergent series, 
what gives 

/(A) = E -7 m ( 1 - ^"(Ao) m^TfiXo)} ■ (3.22) 



„ n\ \ 21 J d\[\ V 21 
In particular 

00 1 / i \ n d n ~ l 
A = Ao + E-777; -JT^ [(^(Ao))1 ■ (3.23) 



{n\ \2tJ dX, 



We accept the formula ( |3.23| ) as the definition of the operator A. Accordingly functions (and 
functionals) of this operator should be understood in the sense (|3.22|) . Obviously, for large value 
of t the difference between Aq and A behaves as t^ 1 — > 0. Therefore we call the functional 
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J-(ip(-),(f),A) asymptotically linear with respect to ip. Nevertheless the replacement A by Ao in 
(|3.19|) is illegal, since ip is unbounded operator. 

Using ( p.23|) one can rewrite the mean value ( |3 . 1 9|) in the form 



(01 £ h ( i)"^r k' (Ao)+ ^ ( ' w ' Ao) ( 1 - ^ I F^w^y-te^o) 



^ n\ \2tJ d)$dp> 
what in turn reduces to the double integral 



13=0 



■ — (0| f°° p dpe 2i{t+i0)p2 I —eP^^'&o+pz^FM-MM+pz) 



n Jo J z 



x (l - ^-V"(Ao + pz)) F c (<f>(v), Ao + pz)e tF «^o+ P z) ^ 



2t 

where the integral over z is taken with respect to a small contour around the origin. Recall that we 
are interesting in the asymptotic behavior of this mean value at t — > oo. It is not difficult to check 
that in this case the term, containing the second derivative tp" , contributes only into a constant 
factor and hence, it can be removed. Then the integrand takes the form considered above, and we 
can average it. We have 



p dpe 2l ^ p2 I —C{p, z, w y F M")M+ P z) ) (3.24) 



-^Tpdpe^U^ 
IX Jo J z 

where 

w(u, p, z) - pz~ l g(\Q + pz, u) - T{i{-,u), w, A + pz) = 0. (3.25) 

It remains now to take the integral in (|3.24j) . For t — > oo it can be estimated via the steepest 
descent method, and we arrive at 

( O | e ^(0,M) Fc ^^^ e tF e (^))| O ) ^ Ce 2it(A-Ao)/3+tF e (^,A,/3),A). (3.26) 

Here the function w(u) satisfies the equation 

w(u, A, /3) - /3£'(A, u) - ^(£(-, u),w, A) = 0, (3.27) 
where numbers A and (3 are defined by the equations of the saddle point: 



2t(3 + -^F e (w,A) =0, 
2*(A-A ) + ^F e (w,A) =0. 



(3.28) 



Thus, in the example considered we again deal with the averaging mapping, however of more 
general form 

<p(v)^w(v), A^A, F(cf),k) -f F(w,A). (3.29) 
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4 Asymptotics of the correlation function 



Now we are ready to consider the asymptotics of the correlation function Ql.lQ . There are two 
different phases in the model of the one-dimensional Bose gas, corresponding to the negative and 
positive chemical potential h in the Hamiltonian ( |1.2|) . While at h < the ground state coincides 
with the bare Fock vacuum, at h > the ground state is the Fermi sphere. At a finite temperature 
the difference is not so significant, but nevertheless the sign of the chemical potential affects on the 
asymptotics of the correlation function. The corresponding formulae were obtained in the paper 
H. Let us present here the explicit expressions. 



;#(o,o)*tOM)) T 



x ex P < 7^— 

Z7T 



(O|C_(0,0|A o )(2t)' iid ^e^( A ) +ltA2 -" A - iM 



2Xt + iip'(X) sign(A - A) In 



_ e <«A)sign(A-A) 

e t + l 



dX |0), (4.1) 



for h < 0, and 



x exp 



2tt 



■ (O|CV(0,0|A o )(2t) 
x - 2Xt + iip'(X) ] sign(A - A) In 



P e ^(Ai)+^(A 2 )+f (A?+Al)-f (Ai+A s )-iW 
e <^>(A)sign(A-A) 



£(A) 

e t 



e r + 1 



-a;(A) 



dX |0), (4.2) 



for /i > 0. 

Let us make necessary explanations. The equations ( [4.1| ), ( [4.2| ) were obtained under the as- 
sumption that t — > oo, x — > oo, while the ratio x/2t = Xq remains fixed. Function e(X) is the 
solution of the Yang- Yang equation (1.4). 

The most important objects for us are three dual fields V'(A), 0(A) and (f>(X). The commutation 
relations between corresponding creation and annihilation operators are 

ic + fi — X\ 



MX),qM] = -MA),5ty(//)] = = In x _ 

[p^(A),^(//)] = \p4>(\),qi>(jJ.)] = M>),qM] = r){\v) = In 



(4.3) 



(A-^) 2 + cV 

As usual all the dual fields commute with each other. Recall that just this property allows one 
at the certain stage of calculations to consider these objects as classical functions. The equations 
( |4.1|) , ( |4.2|) were obtained in H just under such treatment of the dual fields, although some operator 
features were taken into account. In particular, a special attention were paid to the stability of 
the asymptotic formulae under the averaging procedure with respect to the auxiliary vacuum, i. 
e. the corrections, depending on the dual fields, would remain small after their averaging. The 
presence in the asymptotic formula? of the value A, which we dealt already in the previous section, 



A = A + -V / (A), 



(4.4) 
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is caused just by this reason. All the expressions, containing the operator A, should be understood 
in the sense of the series (|3.22|) , (|3.23|) , in particular, the integrals, depending on the sign-function, 
are equal by definition 



/(A,sig n (A-A))dA 



f(X,l)dX 



f(X,-l)dX. 



(4.5) 



One can use ( |3.22| ) for each of the integrals in the r.h.s. of ( f4.5| ). In fact one can say that A arises 
as the result of partial summation of the asymptotic series for the correlation function. 

For operators Ai and A 2 the situation is quite analogous. If we consider dual fields as classical 
functions, then Aj are the roots of the equation 



e{K) =Tsign(A- A i )0(A i 



Ai -> -q T , A 2 



Qt, 



at 



0. 



(4.6) 



where ±qr — zeros of the Yang- Yang function: e(±qT) = 0. In other words the numbers Aj 
correspond to the zeros of the expression, standing under the symbol of the logarithm in the 
integral (O). In the real situation, when 0(A) is the dual field, Aj become operators, defining by 
the series similar to Q3.23j ). We shall simply demand the operator equation ( |4 . 6| ) to transform to 
the usual equation under the averaging mapping 



Ai 



-Qt, 



A, 



Qt, 



at 



0. 



e(Ai) = Tsign(A - Aj)w(Aj), 
Function u(X), entering (|4.2|) , is equal 

ou(X) = sign(A - Ai) sign(A - A 2 ). 



These sign-functions are defined in the same manner as the sign-functions of A ( |4.5| ). 
The equations (|4.1| ), (|4.2j ) contain also a function i^(A), defining the power law of t 



Z7T 



(4.7) 




(4.8) 



Finally the explicit form of factors C±((f), 0|Aq) is unknown (see 0). It is known, however, that 
they do not depend on the field ip, an d they depend on the distance x and time t only through 
the ratio x/2t = Xq, which, recall, remains fixed. 

In spite of the asymptotic formula? ( f4.1| ), ( ^4.2| ) look rather complicated, nevertheless they 
belong to the class of functionals of the dual fields, considered in the previous section. They 
are asymptotically linear with respect to the field ip(X). The operators Aj are some (implicit) 
functionals of the field 0(A). The presence of the third field 0(A), as well as existence of unknown 
factors C±, should not worry us, since due to the trivial commutation relations 
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one can easily see, that the averaging of the field <fi may contribute only into the common constant 
factor. Since our goal is to compute the correlation radius, this factor is not essential, and we even 
can put = 0. 

Thus, applying the averaging mapping (|3.29|) to the equations (|4.1|) , ( |4.2|) and neglecting the 
common constant factor and power law dependency on t, we have 

(*(0,0)*t(z,t)) T — e^*, (4.9) 

where r± correspond to the positive and negative chemical potential respectively. We shall call r± 
correlation radii, although, strictly speaking, they are not necessary real. 
For h< 

_ gtu(A) sign(A-A) 



-t/r_ = 2it(A-X )[3 + itA-ixA-iht + 



2tt 



(x-2Xt) sign(A-A)ln 



e t + l 



dX. 



(4.10) 



Here function w(u), depending on parameters A and /3, is the solution of the integral equation 

g^r — q w W sign (A— A) 



1 f°° 

w(u) = f (A, u) + i{3K(A, u)-— K(u, A) sign(A - A) In 

2,71 J-oo 

where 

£(A,/i) = In 



e r +1 



+ A — /i" 



K{x >» ) = - i ikt {x > ti) = (\-l)* + <*- 



ic + /i — A J 

The parameters A and f3 in turn are defined by two additional equations 



dX, (4.11) 



(4.12) 



i(A - A ) + 



i d r°° 
2^d[3 

1 d 



(A - A) sign(A - A) In 



£{A) 

e t 



w(X) sign(A— A) 



e(A) 

e t +1 



dA = 0, 



+ A _ A ) + — — £jA - A) sign(A - A) In 



e ^p- _ e w(A) sign(A-A) 



e(A) 

e t +1 



rfA = 0. 



(4.13) 
(4.14) 



We also took into account that x = 2tXo. 

For positive chemical potential formula? are rather similar, but there is a set of differences. The 
correlation radius is given by 

it IT 

-t/r+ = 2i*(A - X )(3 +-(A\ + Al- 2h) - - {A x + A 2 ) 



+ 



2tt 



(x - 2 At) sign(A - A) In 



^ _ p w(X) sign(A-A) 



e r 



e t + 1 



-w(A) 



dA. 



(4.15) 



Hereby the function w(-u) satisfies the integral equation 
w{u) = i(£(A ls u ) + £(A 2 , «)) + i/3K(A, it) 



/°° K(u, A) sign(A — A) In 

Z7T J-oo 



_ e «j(A)sign(A-A) 

e t + 1 



dA, 



(4.16) 
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where Ai and A 2 are the roots of the equation fl4.7| ) 



e(A0 =Tsign(A-A i ) M ;(A l ) 



Ax 



A, 



Qt, 



at 



T -> 0, 



(4.17) 



and 

to>(A) = sign(A - Ai) sign(A - A 2 ). 
The parameters A and (3 can be found from additional equations 



1 d r°° 

i{A - A + — ^- / A - A) sign(A - A) In 

Z7T OP J-oo 



i/3 



2ndp 



j_d_ 



e ^T~ — e ^(A)sign(A-A) 

e T + 1 



w(A) 



dA = 0, 



(A - A) sign(A - A) In 



e t 



3 ik(A) sign(A— A) 



e(A) 

e r + l 



-o;(A) 



rfA = 0. 



(4.18) 



(4.19) 



Thus, the exponential decay of the correlation function (1.1) is defined by the function w{u). This 
function can be found from the system of integral equations, which structure is rather similar to 
the Yang- Yang equation ( |1.4| ). Therefore it is quite possible that w(u) can be expressed in terms 
of observables of the one-dimensional Bose gas. In the next section we shall demonstrate this for 
some limiting cases. 

In the conclusion of this section we consider free fermionic limit: c — > oo. Then K(X,fi) = 
£(A, n) = (see (|4.12|) ), and hence, w(u) = 0. Besides e(X) = A 2 — h, A = A , (3 = and 
Ai )2 = ^fy/h. The integral term in the equations ( f4.10|) , ( |4.15|) turns into 



1 

2^ J- 



\x 



2 At I In 



e t 



- 1 



£20 
e t + 1 



dX. 



(4.20) 



In the case of positive chemical potential this expression completely describes the asymptotic 
behavior of the correlation function, since A 2 + A 2 , — 2h = 0. For h < we have additional 
oscillating term —U(Xq + h). These formulae exactly reproduce the results of ||12|| . 

However the free fermionic limit is not too representative. Indeed, in this case all the commu- 
tation relations between operators p and q ( |4.3| ) become trivial, therefore we could put from the 
very beginning all the dual fields equal to zero in the equations (|4.1| ) 



5 Limiting cases 

In this section we consider some limiting cases of the correlation function ( |1 . 1|) . The first case 
corresponds to the autocorrelation. In spite of we considered the limit x — > oo, t — > oo, in fact all 
the asymptotic analysis performed in |Q was based on the assumption that the x/2t = Ao is finite. 
In particular A may be equal to zero, what describes the case x = 0. 

Consider, for instance, the case h < 0. It is easy to check that for Ao = one can put 
A = j3 = in the equations ( 4.11| )-( |4.14j ). Indeed, it follows from ( [4.11 ) that for zero values of 
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these parameters w(u) is an odd function. Differentiating ( |4. 1 1| ) with respect to A and (3 at zero 
point we find that the derivatives dw/dA and dw/dj3 are even functions of u. Then the equations 
(|4.13|) , (|4.14j) are valid automatically, and the result slightly simplifies 



ih 



2 

7T JO 



A In 



e -^r L _ e w(\) 

e t + l 



dX, 



(5.1) 



where w(u) is the solution of the equation 

1 



w[u) 



2tt Jo 



(if(A, u) — K(X, —u)J In 



e(A) 

e t 



e t + 1 



dA. 



(5.2) 



Similarly for positive chemical potential we obtain 



r; 1 = - A^) 



7T JO 



A In 



£P0 

e t 



e r + 1 



-a;(A) 



G?A, 



(5.3) 



and 



«;(u) = -(e(A 1 ,«)-e(Ai,-«)) + ^ 

2 27T Jo 



K(A, u) - K{X, -u)) In 



e A F L _ e w(\) 

e t + 1 



w(A) 



dA. (5.4) 



Hereby it turns out that A 2 = —Ax. 

The most interesting is the low temperature limit for the case h > 0. In this limit we have a 
possibility to compare our results with the known ones. Indeed, as we have mentioned already, 
at T = the asymptotics of the correlation function must be power-like. This means that the 
correlation radius r + — > oo. Besides, the correlation radius was computed for small temperature 
in the linear approximation in M. Thus, one can consider the low temperature limit as a good 
test for our results. Besides we have a possibility to make more precise estimates of the results 
obtained previously. 

First, let us present the system of equations, describing the ground state for h > [pL3[| . Recall 
the Yang- Yang equation (|1.4j) 

T r°° / 
e(X) = X 2 -h / K(A,/i)ln 1 + e' 

277 J-oo V 



£(M) 
T 



dfi. 



(5.5) 



The function e(X) has two real roots e(±g r ) = 0, and e(A) > 0, if |A| > q?, and e(A) < if 
|A| < qr- At T — > the roots ±qr tend to some fixed values (corresponding to the boundaries of 
the Fermi sphere): q? — > q- It is easy to see that 



exp 



T 



T^O 



0. 

CXO, 



\H\ > q, 
< q. 
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Thus the equation (513) turns into 

e (A) = A 2 - h + — f K(X, fi)e (n) d/x, e (±q) = 0. 

Z7T J-o 



(5.6) 



Here and further we denote observables at zero temperature by the subscript 0. Similarly the 
equations for the density and the scattering phase have the form 



Po(A) 



K(X, /i)po(/i) dfi, 



2tt 2tt J- q 

F {\, u) = -L f (A, /i)F (/i, i/) d/x + -^£(A, i/) 
The derivative of the momentum and the velocity as before are equal to 

dk Q (X) £n(A) 

-27T Po (A), u (A)- 



<9A 



(5.7) 
(5.8) 

(5.9) 



Let us give also the equation for the resolvent of the integral operator, entering the equations 
O-fJ): (I + R)(I-±K) = I 



R(X, v) = ^-l q K(X, ijl)R(ji, u) d(i + ^-K(X, v). 

271 J-q 271 



(5.10) 



It is easy to see that R(X, /i) = —d fl F (X, p,). Besides, comparing the equations ( |5.7| ) and ( |5.8| ) and 
using K(X,n) = —id\£(X, n), we find: 



27Tp (A) = 1 + F (X, -q) - F (A, q). 



(5.H) 



Similarly to the Yang-Yang equation the integral equation for the function w(u) becomes linear 
in the limit T —> 0. Let, for example, Ai < A2 < A. Then instead of ( 4.16 ) we have at T = 

Mu) = Uti-q, u) + t{q, u)) + i/3K(u, A) + -L f K(u, X)w (X) dX. (5.12) 

I 2lX J-q 



Here we have used that A2 = — Ai = q for T = 0. Comparing the equation ( 5.12 ) with (|5.8|) and 
(15T0D , we find 

w (u) = -i7r(F (u,-g) + F (u,g)) + 27ri^i2(u,A). (5.13) 
The restrictions for the numbers A and (3 have the following form: 

1 d 



i(A — Aq) — 



2tt<9/3 



(A - A)^o(A) dA = 0, 



i/3 



J_d_ 

2^dA 



(A - X)wq(X) dX = 0. 



(5.14) 
(5.15) 
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It follows immediately from ( p-13| ) and ( 5.15| ) that (3 = 0, and hence, 

w (u) = -iir(F (u, -q) + F (u, q)) (5.16) 
Substituting this into ( 4.15[) (where the logarithm also turns into linear function), we find 

- t/r+ = tt(q 2 ~h) + % - J" (x - 2Xt)(F (X, q) + F (A, -q)) dX. (5.17) 



It is easy to see that the r.h.s. of (|5.17|) is equal to zero identically. Indeed, due to ( |5.8| ) the 
integral in the r.h.s. of (|5.17 ) can be represented in the form 

— J\x - 2Xt){5{X -fi) + R{X, //))(£(//, q) + £(//, -?)) dXdfi. 

Acting by the resolvent to the left and using (|5.6|) - (|5T9|) , we obtain 

i J\xk' (X) - te' (X))m, q) + e(A, -q)) dX. 

Since k' (X) is an even function, then the coefficient at x is equal to zero. The remaining integral 
after the integration by parts gives 

£; J*e (X)(K(X,q) + K(X,-q))dX. 

Now one can use the equation (|5.6|) and the condition Eo(±q) = 0, after that we obtain 

for T = 0. (5.18) 

Thus our formulae do give the correct result at zero temperature. 

Recall that the introducing of the dual fields is equivalent in some sense to factorization of 
the S'-inatrix. It is natural to expect that under the averaging of the dual fields the objects have 
to appear, having direct relation to the scattering matrix. However, while the factorization takes 
place at the bare Fock vacuum, the averaging is performed already in the thermodynamic limit, 
i.e. at the physical vacuum. It is remarkable that at this rate instead of the 'bare' S'-matrix (|2.3|) 
the function wq(u) arises, which can be expressed in terms of the 'dressed' scattering phase (|5.8|) 
at the boundary of the Fermi sphere. 

Observe also, that in the calculations presented above, we did not need the explicit value of A. 
However, it is not difficult to find it. It follows from the equation ( p-14| ) that 



A -A c 

and we immediately find 



f q (X - X)R(X,A) dX = 0, 

J— a 



Xok' (A) - y o (A) = 0, or «b(A) = |. (5.19) 
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Thus the quantity A has a pure physical sense: a particle possessing velocity v (A), passes the 
distance x spending the time t. 

In order to compute the correlation radius for small, but not zero temperature, more accurate 
estimates of the integrals containing logarithmic function are necessary. The example of such 
calculations is given in the Appendix A. It is clear, however, that all low temperature corrections 
to the correlation radius, as well as to the integral equation for the function w(u), are provided 
by vicinities of the points ±g, defining the boundary of the Fermi sphere. Let us present here the 
integral equation for the function w(u), arising in the first order with respect to T. We shall look 
for the function w(u) as a series 

w(u) = w (u) + Twi(u) + T 2 w 2 (u) + . . . 

Here Wo(u) is given by ( |5.16| ). For Wi(u) we have 

Wl (u) - -L J* K(u, AWA) dX = -E±^( Wo (q) + ,vr) 2 - (w {-q) - ztt) 2 . (5.20) 



Hence 



7T 2 



wi(u) = ^773(1 - Fo(Q, Q) - Fo(q, -q)T(R(u, q) + R(u, -q)). (5.21) 

Here we took into account, first, that F (— A, — /i) = — F (X,fi). Second, we have put already the 
parameter (3 = 0. On can easily check that it is always so in the framework of the low temperature 
approximation. Indeed, for any order of T the dependency of w(u) on the sign-functions sign(A— u) 
disappears. The parameters A and j3 enter the equation for w(u) only in the combination (3K(A, u) 
(see ( |A.8| )). Therefore the derivative dw/dA always turns out to be proportional to (3, and hence, 
the equation (|4. 1 9|) is valid for (3 = 0. 

Substituting w\{u) into the answer for the correlation radius, after simple algebra we arrive at 

r; 1 = — (1 - F (q, q) - F (q, -q)f{\x - v t\ + \x + v t\). (5.22) 

Here Vq = v (q). This result is given for arbitrary location of the roots A x and A 2 with respect 
to the point A. In the case Ai < A2 < A, which we considered from the very beginning, one can 
remove the symbols of the absolute value. 

Comparing ( |5.22|) with the result obtained in 

_i _ *T \x-v t\ + \x + vpt\ 

+ ~ 4^o ' (27Tp (g)) 2 ' { } 

we obtain an identity 

-\- = 1 - F {q, q) - F (q, -q). (5.24) 
27rp (g) 

At the first sight such the non-linear relationship between the solutions of the linear integral 
equations looks rather strange. However, taking into account a non-linear identity for the scattering 
phase 0, [[L5J 

(l-Fo(g,g)) 2 -F 2 (g,-g) = l (5.25) 
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and the equation ( 5.11j ), we make sure of the validity of the identity ( |5.24j ). Thus, in the linear ap- 
proximation with respect to the temperature our result for the correlation radius exactly coincides 
with the result obtained in the framework of the conformal field theory. 

The computation of higher temperature corrections does not meet serious difficulties. In the 
all orders T n the equations for w n (u) remain linear (see ( |A.8| )), and their solutions can be written 
explicitly in terms of the resolvent of the operator / — ^-K and its derivatives. Respectively the 
correlation radius is expressed in terms of functions k' Q (q), e' (q), Fo(q,±q) and their derivatives. 
In particular the correction to r^ 1 of the order T 2 vanishes. Thus the result ( |5.22| ) is valid up to 
the terms of order T 3 . 



Conclusion 

We have considered the large time and long distance asymptotic behavior of the correlation function 
of the quantum one-dimensional Bose gas. As we have shown, the asymptotics can be expressed 
in terms of the integral equation solutions, which are closely related to the equations of the 
thermodynamic Bethe Ansatz. In the low temperature limit we were able even to describe the 
solutions of mentioned equations in terms of observables of the model. We would like to draw 
attention of the reader to the fact, that these integral equations arise under the averaging of the 
dual fields. While in the earlier works these quantum operators played purely auxiliary role and 
mostly were used for solving of certain combinatorial problems, the results obtained now allows 
us to have a look at these objects from a new point of view. In particular, it is hardly accidentally 
that at the zero temperature the solutions of the integral equations depend on the scattering phase. 
Apparently the role of the dual fields is much more fundamental, than it was understood before. 
We would like to thank A. R. Its and V. E. Korepin for useful discussions. 



A Appendix 

As an example of the low temperature expansion we consider the integral, entering the Yang- Yang 
equation 

I = T 1°° K{X,fi) In (l + e-^ dfi. (A.l) 
Obviously, at T = only the interval [— q, q] contributes into this integral 

I T=0 = - [ K(X,n)e(ji) dfi, e(±q) = T = 0. (A.2) 

J-q 

At small, but finite temperature the corrections to this expression arise due to the contributions 
of vicinities of the points ±q. Consider, for instance, the integral 

/ + = Ty ff(A, //) In (l + e~T- J d/i, (A.3) 



20 



where e{qx) = 0. Making the replacement of variables e(fi) = Tz, we obtain 

I + = T 2 j™ (p(X, Tz) In (l + e~ z ) dz, (A.4) 

where tp(X, Tz) = K(X,fi)/e'(ii). The integral flPD can be expanded now into the series with 
respect to T. We have 

oo 

I + = £ T" + V (n) (A, 0)(1 - 2-«- 1 )C(n + 2), (A.5) 

n=0 

or in terms of the functions K(X,fi) and 

) C („ + 2 Wt ^M, where D x = ^. (A.6) 

Similarly all the remaining contributions from vicinities of the points ±qx can be computed. Then 
one need to substitute the expansion obtained into the Yang- Yang equation and to put 



e(X) = J2T n e n (\), 

n=0 

(A.7) 

oo 

q T = q+Y. Tn( ln- 



n=l 



The arising infinite set of equations can be solved via recursion. As a result the functions e n and 
parameters q n are expressed in terms of e' (q) and the resolvent R(X, ±q). 

Similarly all the integrals of the section 4 can be decomposed into series. Let us present the 
complete expansion with respect to the temperature of the integral equation ( [4.16D for the function 

w(u) 

w{u) - ±- [ q K(u, X)w(X) dX = \ (a-q, u) + £(q, u)) + if3K(A, u) 

2tt J-q 2 v ' 



J_ « (-2^ 7 )"T"- 1 Dn _ 2 K(u, igT 

**4=kh* n] qT z'M 



B n (l +1 ^)-B n{ l/2) 



. (A.8) 



Here B n (z) are the Bernulli polynomials. It is yet necessary to substitute into this equation 
decompositions ( |A.7|) , after that functions w n {u) can be evaluated in terms of s' (q), R(X, ±q) and 



their derivatives. 
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